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0\ ■ Abstract 

' The properties of relativistic particles in the quasiclassical region are investi- 

gated. The relativistic semiclassical wave equation appropriate in the quasiclassical 
^ , region is derived. It is shown that the leading-order WKB quantization rule is the 

I appropriate method to solve the equation obtained. 

00 ; PACS number(s): 03.65.Ge, 03.65.Sq 

1. Introduction 
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It is generally agreed that the full description of quantum systems is given by quantum 
O ' mechanics. However, investigations by many authors show that, in fact, there exists an 

alternative approach. This is the approach based on the quasiclassical method. 
Q\ . The general form of the semiclassical description of quantum-mechanical systems has 

^ ! been considered in Ref. |Q . It was shown that the semiclassical description resulting from 

Weyl's association of operators to functions is identical with the quantum description 
and no information need to be lost in going from one to the another. The semiclassical 
^ ■ approach merely becomes a different representation of the same algebra as that of the 

quantum mechanical system, and then the expectation values, dispersions, and dynamics 
Q^! of both become identical. What is more "the semiclassical description is more general 

than quantum mechanical description..." ||l|]. 

The quasiclassical method is based on the correspondence principle between classical 



' functions and operators of quantum mechanics and with a real, but not necessary posi- 



es 



tive, probability density function in phase space corresponding to a particular quantum- 
mechanical state. The correspondence principle is used to derive the wave equation in 
quantum mechanics. In Ref. this principle has been used to derive the semiclassical 
wave equation appropriate in the quasiclassical region. It was shown that the standard 
WKB method (to leading order in h) is the appropriate method to solve this equation. 
Solution of the obtained wave equation by this method yields the exact eigenvalues for 
all solvable spherically symmetric potentials. The corresponding eigenfunctions have the 
same behavior as the asymptotes of the exact solutions. 

In this work we develop the semiclassical approach for relativistic systems. We begin 
our analysis with the relativistic classic consideration of the problem in the Hamilton- 
Jacobi formulation. Two different versions are considered: in the first the potential is 
zero component of a four-vector, and in the second the potential is Lorentz-scalar. In 
the second version, the total energy of the system is presented by analogy with non- 
relativistic case as sum of the kinetic energy and the potential energy of interaction. 
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Resulting equation is written in covariant form of the relativistic one-particle equation. 
After that, using the correspondence principle, we derive a relativistic semiclassical wave 
equation appropriate in the quasiclassical region. 

It is shown that the leading-order WKB quantization condition is the appropriate 
method to solve this equation for the spherically symmetric potentials. Solution of this 
equation for the Coulomb and linear potential is given. In particular, we show that unlike 
the existing relativistic wave equations, solution of the semiclassical wave equation for the 
scalar-like potential results in the eigenfunctions which are regular at the spatial origin. 

2. Relativistic Semiclassical Wave Equation 

There are several wave equations to describe the motion of relativistic particles: the 
Klein-Gordon equation, the Dirac equation, the Bethe-Salpeter equation, the spinless 
Salpeter equation, quasipotential equations. These equations can be solved exactly for 
some spherically symmetric potentials. However, S'-wave solutions of these equations, 
for example, for the Coulomb potential diverges at the spatial origin. Wave-function 
divergence at the spatial origin is actually a general problem affecting relativistic wave 
equations. For example, the solution of the Dirac equation for the Coulomb potential for 
the S'-wave states behaves as oc t-v^i-^^-i at r and the divergence of the solution 
to the spinless Salpeter equation is much stronger: ip oc r"^""/^'^ 

Essential problem of relativistic description of physical systems is concerned with 
the entity of the potential. The potential is a non-relativist ic conception and it has been 
correctly defined in classic mechanics as the function of the spatial variables. In relativistic 
case the potential has not been correctly defined. A serious problem is the nature of the 
potential: whether it is Lorentz-vector or Lorentz-scalar or it is their mixture? This 
problem is especially important in hadron physics where, for the vector-like confining 
potential, there are no normalizable solutions In this letter we investigate both of 
these possibilities when the potential is the Lorentz-scalar and zero component of a 4- 
momentum. 

To derive the wave equation with the help of the correspondence principle, we need to 
know the corresponding classic equation. This is why we consider first classic description 
of relativistic particles. The classic equations are then used to derive the corresponding 
semiclassical wave equations. 

As known, for a relativistic particle moving in the field of the spherically symmetric 
potential V{r), the interaction is introduced according to the gauge invariant principle, 

^ P^= p^ — eA^, where is the four-momentum of the particle and is the external 
field. For zero component Pq = E this results in the total energy of the particle in the 
form (hereafter we basically use the system h = c = 1) 

E = + + V{r), (1) 

where V{r) = cAq is the static potential. Equation (1) (with the help of the correspon- 
dence principle) results in the spinless Salpeter equation |]5|. For the classic equation 
p'^ = [E — V^(r)]^ — m^, the same correspondence principle results in the static Klein- 
Gordon equation. 
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Due to the success of the non-relativistic potential models, it can be used to scrutinize 
the non-relativistic method. This method has proven extremely successful for the descrip- 
tion not only of non-relativistic systems but also of relativistic bound states. Potential 
models work much better than one would naively expect (see Ref. 0). This success 
is somewhat puzzling in that it persists even when the model is applied to relativistic 
systems. 

It is generally agreed, that in nonrelativistic classic mechanics, the total energy E' of 
a particle of mass m moving in the field of the central-field potential V{r) is given by 
equation E' = T + V{r), where T is the kinetic energy of the particle; the total energy E' 
is the integral of motion. In relativistic case, the total energy can be written in analogous 
form and, what is more, the potential V{r) in this case is a Lorentz-scalar. 

Energy a/pM-^ttS' Qf relativistic particle in Eq. (1) can be presented as sum of the 
kinetic energy T and the rest mass m as y/p"^ + m? = T + m. Then equation (1) takes 
the form 



E = T + W{r), (2) 

where we have introduced notation W{r) = m + V{r). From nonrelativistic point of view, 
W{r) is the potential energy shifted by constant value m with respect to V{r). On the 
other hand, from relativistic point of view, W{r) can be considered as a variable mass 
(this topic has been discussed in series of works 0, H, ^). This makes sense because two 
ingredients of the total energy (2), the rest mass m and the potential V{r), are static 
quantities and they can be combined into one object, W{r). Thus we can write for the 
kinetic energy T of the particle of the effective mass W{r), T = E—W{r) = tt"^ /[E+W{r)], 
where the squared momentum vf^ is 

n^ = E^- W\r). (3) 

Equation (3) has the form of the equation for a free relativistic particle with the effective 
mass W{r) and can be written in the covariant form as 

ttX- W^'W = 0, (4) 

where n'^ = {E, vf). 

Now, we can derive the wave equations which correspond to the classic equations 
written above. We start with the classic problem in the Hamilton- Jacobi formulation. 
The static Hamilton- Jacobi equation corresponding to the classic equation p"^ = [E — 
V{r)]'^ — w? is 

where 5*0 is the action of the system. To obtain the corresponding wave equation we 
present the wave function in the form ipif) = A exp ^S{r) , where A is the arbitrary 
constant. As was shown in Ref. |^, in the quasiclassical region, in representation of the 
wave function '?/'(r*), the squared momentum operator has the form 
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P' 



-ih 



dr 



+ 



-ih 



de 



+ 



sin^ 6 



-ih 



(6) 



which follows from the classic expression for with the help of the correspondence 
— — ^ — I'ft-^ " — r,9,Lp p|. Thus we obtain the wave equation {h = c = 1) 



principle ^ 




^{r)= {E -V{r)f -m^ i){f). (7) 



Equation (7) is the second-order differential equation in canonical form. This relativistic 
semiclassical wave equation is closely related to the classic equation (5) and is appropriate 
to describe relativistic quantum systems in the quasiclassical region. Below, we solve this 
equation for some spherically symmetric potentials. 

Let us consider now the classic equation (3). The corresponding semiclassical equation 
can be derived analogously to the above case. Using the operator (6), with the help of 
the correspondence principle we obtain the following equation 



dr 



+ 



''de 



1 



sin 9 



ip{r). 



The correlation of the function il){r) with the wave function ip{r) in case of the spher- 
ical coordinates is given by relation tp{'r) = det gijip (r ) , which follows from the iden- 
tity: / I ip{r) p d^r = J I ilj{r) ^ det Qijdr d6 dip = 1, where Qij is the metric tensor 
(det (yfjj = sin 9 for the spherical coordinates). 



3. Solution of the Semiclassical Wave Equation 



Relativistic wave equations are usually solved in terms of special functions, with the 
help of specially developed methods or numerically. However, almost together with quan- 
tum mechanics, the appropriate method to solve the wave equation has been developed; 
it is general simple for all the problems, and its correct application results in the exact 
energy eigenvalues for all solvable potentials. This is the phase-integral method which is 
also known as the WKB method |Tl| . 

The WKB method in its present form has been formulated to solve one-dimensional 
two-turning-point problems. Within the framework of the WKB method the solvable 
potentials mean those potentials for which the eigenvalue problem has two turning points. 
However the WKB method can be used to solve problems with more then two turning 
points. 

The WKB formulas can be written differently, i.e. on the real axis [l^ and in the 
complex plane [|1T|. For the two-turning-point problems, the standard leading-order WKB 



quantization condition is defined by 



^2 



p^[x)dx = Tihyn + ^ 



(9) 
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where p{x) is the generahzed momentum and n is the number of zeros of the wave function 
inside the interval [xi,a;2]. However, the most general form of the WKB solution and the 
quantization condition can be written in the complex plane. 

The WKB solution of the one- dimensional Schrodinger's equation is searching for in the 
form ip{z) = aexp -^S^z) , where a is the arbitrary constant and the action of the system, 

S{z), is written as the expansion in powers of h, S{z) = So{z) + hSi{z) + h'^S2{z) + . . .. 



In the leading h approximation we have Sq{z) = ± p{z)dz, Si{z) = ifiln ^J\p{z)\, where 
p{z) is the classic momentum. Thus, in this approximation, the WKB solution is 



— (^) = ^_exp 

'\p{z)\ 



lfp{z)dz_ 



(10) 



In quantum mechanics, quantum numbers are determined as number of zeros of the 
wave function in the classically allowed region. In the complex plane, number of zeros 
of a function y{z) inside the contour C is defined by the argument's principle [Q. For 
the wave function '^/'(-z), according to this principle we have 

/ = 2mN, (11) 

Jc ipiz) 

where ip/{z) is the derivative of the function ip{z). Contour C is chosen such that it 
includes cuts between the turning points on the generalized momentum p{z) > and no 
other singularities. 

Substituting (10) into (11) we obtain the quantization condition 

^ ^p^{z)dz + i^j> ^^dz = 2'KnN. (12) 

In the case, when p{z) is a smooth function of the spatial variable and equation p'^{z) = 
has two roots (turning points), Eq. (12) results in the quantization condition 

j^p'^{z)dz = 27ch (^+ ^) • (13) 
In particular case p'^{z) = const quantization condition has the form 



\ljp\z)dz 



2nhN. (14) 



Quasiclassical method is fruitfully used not only in quantum mechanics but also in the 
relativistic theory [jl3|, [1^, 0. Application of the leading-order WKB quantization rule 



to the nonrelativistic semiclassical wave equation p[ yields the exact energy eigenvalues 
for all solvable spherically symmetric potentials and no further Langer-like corrections are 
necessary. 

Equations (7) and (8) are the second-order differential equations of the Schrodinger 
type in canonical form. Important feature of these equations is that, for the two and more 
turning-point problems, they can be solved exactly by the conventional WKB method. 

For the spherically symmetric potentials, Eqs. (7), (8) are separated. Angular equation 
for each of them is the same, 
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d 



1 



d 



89 j sin^ 9 \ dip ^ 
and the corresponding radial equations are 



Y{9,^) = M^Y{9,^), 



(15) 



—I 



d 


)■ 


dr 




—i 


d 


dr 



{E-V{r)f-m^--^ 



3.1. The angular momentum eigenvalues 



M2 



(16) 
(17) 



Angular equation (15) determines the squared angular momentum eigenvalues 
which enter into the radial equations (16), (17). The solution of equation (15) has been 
obtained in Ref. B that gives for M^, 



/ = 0,1,2, 



This same result can be obtained in the complex plane with the help of the quantization 
condition (13). It is interesting from the methodological point of view to bring here the 
calculation of in the complex plane because this calculation is shorter and beautiful. 
Separation of the angular Eq. (15) results in two one-dimensional equations. 



-4) 



dip 



Ml 

sin^ 9 I 



(19) 
(20) 



The WKB quantization condition (14) which is appropriate to the angular equation (20) 
[p(y9) = Mz], gives Mz = m, m = 0,1,2, .... The corresponding quasiclassical solution is 



m(V5) 



■imip 



(21) 



where A and B are the arbitrary constants. 

The WKB quantization condition (13) appropriate to equation (19) is 



/ - M2 
'M2 ^d9 



27T [ne + 



= 0,1,2, 



(22) 



ic V sin^^ " V" ' 2. 

To calculate integral (22) (as other hereafter) we use the method of stereographic projec- 
tion. This means that, instead of integration about a contour C enclosing the classical 
turning points, we exclude the singularities outside the contour C, i.e., at 6* = and oo in 
this particular case. Excluding these infinities we have, for the integral (22), I = Iq + I^. 
Integral Jq = —2ttMz, and I^o is calculated with the help of the replacement z = e*^ that 

gives loo = 2n\f^^. Therefore / = 2n{\/ — M^) and we obtain eigenvalues (18). 
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Equation (18) represents the squared angular momentum eigenvalues in the quasiclas- 
sical region which are different from the known ones l{l + 1) obtained from the solution 
of the Klein-Gordon equation. The centrifugal term in the radial equations (16), (17) has 
the form (/ + |)^/r^ for all spherically symmetric potentials V{r) which is not equal to 
zero at / = 0. The WKB solution corresponding to the eigenvalues (18) has the correct 
asymptotic behavior at 6' — > and tt for all values of / and agrees with the asymptote of 
the spherical functions Yim{0, 0- It is important to emphasize that the squared angu- 
lar momentum eigenvalues (18) have obtained in our approach from the solution of the 
angular semiclassical equation (15) in the framework of the same quasiclassical method. 



3.2. The Coulomb potential V{r) = —a/r 



Consider several examples to show the efficiency of the semiclassical approach in rel- 
ativistic theory. Let us deal with first the classic problem for the Coulomb potential, i.e. 
the radial semiclassical equation (16). The quantization condition (13) appropriate to Eq. 
(16) with the Coulomb potential is 



, 2aE A2 / 1\ 

. E^-m^ + - = 2n(nr + -), 23 

c V ^ V 2/ 

where = (/ + 1)^ — and a contour C encloses the classical turning points ri and r2. 
Using the method of stereographic projection, we should exclude the singularities outside 
the contour C, i.e. at r = and oo. Excluding these infinities we have, for the integral 
(23), / = 2n{aE/V-E^ + - A), and for the energy eigenvalues this gives 

m , , 

Er, = , =. 24 

- _| a'' 

(n.+i+A)' 

Equation (24) coincides with the one obtained from the solution of the Klein-Gordon 
equation with the Coulomb potential by known exact method. 

Now, let us deal with the semiclassical wave equation (17). The quantization condition 
(13) appropriate to Eq. (17) with the scalar-like Coulomb potential is 



The phase-space integral (25) is calculated analogously to the above case in closed form 
that results in the energy eigenvalues 



Er,. = m 



+ 1. (26) 



Formula (26) can be written in apparent relativistic form, E\ = m? + p^, where p„ is 
given by equation 

Pn = , (27) 

+ i + J{1 + i)2 + a2 
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The radial quasiclassical eigenfunctions inside the region of classic motion [ri,r2] are 
written with the help of the Eq. (10) 



i?^^«(r) = -^== cos (^f^ ^^)dr - . (28) 



'\p{.r)\ 

So far, as \p{r)\ ~ - at r — 0, this gives, for the WKB solution in the representation 



of the wave function ■?/'(r), R{r) = R^^^{r)/r oc r^/^. Because A = \J {I + \Y + ^fl) 
for the S-wave states, this results in R{r) oc r", i.e. R{r) is regular at small r. We 
obtain this result because the potential in Eq. (17) is scalar-like. Note that in the case 
of the Klein-Gordon equation the radial S'-wave function R{r) diverges at r — 0, i.e. 
R{r) oc r~". 

The quasiclassical solution (28) can be written in the form of a standing wave. For 
this note that the hermiticity of the operator (6) implies the adiabatically slow alteration 
of the derivatives q = r,9, Lp, i.e. ^ ~ const 0]. This condition anticipates the final 
result, i.e. discrete constant eigenvalues p„ of the operator (6). Integration of equation 
^ ~ const gives, for the action 5*0, So{r) = \pn\r + const, where p„ is given by Eq. (27) 
(for this problem). In the region of classic motion, where p{r) > 0, the solution has the 
form of a standing wave 

R^{r) = C cos (\pn\r-xi -I), (29) 

where xi is the value of the phase integral (25) at the turning point ri. Analogous 
solutions can be written for other one-dimensional equations obtained after separating of 
Eq. (8). 

Thus we obtain two results for the Coulomb potential by the WKB method: the known 
exact eigenvalues for spinless particles (24) which coincide with those obtained from the 
solution of the Klein-Gordon equation and the other ones (26) obtained from the solution 
of the semiclassical equation (8) with the scalar-like Coulomb potential. We compare two 
relativistic formulas (24) and (26) for the energy eigenvalues of an electron moving in the 
field of the static spherically symmetric Coulomb potential in Table 1. 
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Table 1 . Energy levels of the hydrogen atom 
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-0.2127435 


-0.2127430 



In Table 1 denotes the eigenenergies obtained from solution of the Schrodinger 

equation and E^'^ are the energy eigenvalues obtained from the Klein-Gordon equation, 
I and rir are the orbit and radial quantum number, respectively. Eigenenergies E^^ have 
been calculated with the help of Eq. (26). 

As is well known, the energy spectrum resulting from the Klein-Gordon equation with 
the Coulomb potential contains a fine-structure effect giving an over-estimation of the 
fine-structure level spreading, while the result from non-relativistic Schrodinger equation 
does not [|1^]. In our case, the eigenvalues E^'~^ are lower than E^^ and E^^ , hereby 
E^^ are closer to non-relativistic ones. Fine-structure term which follows from E^^ is 

where n = rir + I + 1. This term is different from the corresponding correction which 
follows from E^'^ by sign and the second term. 



3.3. The linear potential V{r) = nr 



The linear potential is one of the special interest in hadron physics. In relativistic 
case, this potential results in the linear Regge trajectories which observed at experiment. 
In relativistic potential models of quarkonia based on a Dirac-type equation with a local 
potential is a sharp distinction between a linear potential V which is vector-like and 
one which is scalar-like. There are normalizable solutions for a scalar-like V but not 
for a vector-like V Q]. No other problems arise and no difficulties encountered with the 
numerical solution if the confinement potential is purely scalar-like. It was shown in many 
works that the effective interaction has to be scalar in order to confine particles inside the 
hadrons (see, for example, Refs. @, |16|)- 

In relativistic theory the linear potential represents the so-called "insoluble" problem 
with four turning points. This problem is insoluble also from the viewpoint of the conven- 
tional WKB method. However, this problem can be solved by the WKB method in the 
complex plane. In the semiclassical equation (17), the potential is a scalar-like. Consider 
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this equation for a two-particle system of equal masses interacting by means of the linear 
potential. The radial semiclassical equation (17) is then: 



2 



4) 



, ^ >2 {l + k 

— [m + Kr 



4 r 



2 



R{r), (31) 



From four turning points of this problem, two of them, ri and r2, lie in nonphysical region 
r < 0, and the other two, and lie in the physical region r > 0. To calculate the 
phase-space integral in the complex plane we chose a contour C enclosing the cuts (and, 
therefore, zeros of the wave function) at r > and r < between the turning points ri, 
r2 and r^, r^, respectively. Thus, the quantization condition appropriate to Eq. (31) is: 

^Tiirir + l), (32) 




/c V 4 ' ' r2 V 2. 

where denotes zeros of the wave function at r > 0. 

Outside the contour C, the problem has two singularities, i.e. at r = and oo. Using 
the method of stereographic projection, we have, for the integral (32), / = Jq + /oo, where 
Jo = —2-n{l + i) and the integral I^o is calculated with the help of the replacement of 
variable, i.e. = 7, that gives I^o = 2ttE'^ / {8k). Therefore, for E"^, we obtain 

= 8k (2nr + 1 + 1), (33) 



2 

which is very similar to that of a harmonic oscillator-type Hamiltonian. 

It is an experimental fact that the dependence E^{1) is linear for light mesons (linear 
Regge trajectories [|l3])- However, at present, the best way to reproduce the experimental 



masses of particles is to rescale the entire spectrum assuming that the masses M of the 
mesons are expressed by the relation |T^ 



M„2 = El- C\ (34) 

where En is given by (33) and C is a constant energy (free parameter). Relation (34) 
is used to shift the spectra and appears as a means to simulate the effects of unknown 
structure approximately. 

However, the required shift of the spectra naturally follows from the solution of the 
semiclassical equation (17). To describe the Regge trajectories observed at experiment 
one needs to take into account the week coupling effect. In this case the potential includes 
the Coulomb-like term, — ag/r, i.e. the scalar-like potential is V{r) = —ag/r + nr (the 
funnel potential). For light mesons one may expect that the qq bound states will feel only 
the linear part of the potential which gives the main contribution to the binding energy. 
We thus assume that the Coulomb term can be considered as a small perturbation. Then 
we obtain instead of (33) 

El = 8k (2nr + l-a, + ^) . (35) 



2. 

Equation (35) does not require any additional free parameter. We obtain the nec- 
essary shift with the help of the term — S/ta^ which is the result of interference of the 
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Coulomb and linear terms of the interquark potential. Note that we obtain the correct 
sign (minus) for this term only in the case of the scalar-like potential. Equation (35) gives 
a good description of mass of light vector mesons for the typical value n ~ O.lAGeV^, 
i.e. reproduces the linear dependence of angular momentum / on E'^ that is in agreement 



with the experimental data for light mesons |17 



4. Conclusion 



In this letter, we have developed the semiclassical approach to the description of rela- 
tivistic systems. In our analysis we have started with the relativistic classic consideration 
of the problem in the Hamilton- Jacobi formulation. Two different possibilities have been 
considered: in the first the potential is the zeroth component of a four-vector and in the 
second, the potential is the Lorentz-scalar. 

The main result is the derivation of the relativistic semiclassical wave equation [(7) 
for the vector-like potential and (8) for the scalar-like potential] appropriate in the quasi- 
classical region. This equation is the second-order differential equation in canonical form 
(without first derivatives). In case of the scalar-like potential, the resulting equation has 
been written in covariant form. Unlike the Klein-Gordon equation, the semiclassical one 
(in case of the vector-like potential) results in another integral of motion, i.e. the squared 
angular momentum = (/ + |)^/i^. This is same squared angular momentum we have 
obtained also for the semiclassical equation with the scalar-like potential. 

The squared angular momentum eigenvalues = {I + |)^^^ have been obtained in 
our approach in a natural way from the solution of the angular semiclassical equation 
(15) in the framework of the same WKB method that gives the justification of the Langer 
correction as that to the squared angular momentum eigenvalues. This means that the 
centrifugal term in the radial semiclassical equation has the form (/ + |)^/r^ for any 
spherically symmetric potential V{r). 

Important feature of the semiclassical wave equation is that it can be solved by the 
WKB method. We have shown that the standard leading-order WKB quantization con- 
dition reproduces the exact eigenvalues for two-turning-point problems. Besides, this 
equation can be solved by the WKB method for more then two turning-point problems. 

To show efficiency of the semiclassical approach in relativistic theory we have con- 
sidered several examples. In particular, we have obtained two results for the Coulomb 
potential by the WKB method: the known exact result for spinless particles (24) which 
coincides with one obtained from solution of the Klein-Gordon equation and another re- 
sult (26) obtained from the solution of the semiclassical equation (8) for the scalar-like 
Coulomb potential. We have shown that, unlike the known relativistic wave equations 
for the Coulomb potential, the semiclassical one for the scalar-like Coulomb potential 
has the regular solution at the spatial origin. Conception of the scalar-like potential is 
especially important in hadron physics. The solution of equation (8) for light qq states 
with the funnel potential has allowed us to reproduce the linear dependence £"^{1) and 
obtain apparent form of the shift parameter C^, = —San. 

Acknowledgment. This work was supported in part by the Belarusian Fund for Fun- 
damental Researches. 
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